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ABSTRACT 

As part of a multi-faceted program to investigate the origin and early evolu- 
tion of sub-stellar objects, we present high-resolution Keck optical spectra of 14 
very low mass sources in the IC 348 young cluster and the Taurus star-forming 
cloud. All of our targets, which span a range of spectral types from M5 to M8, 
exhibit moderate to very strong Ha emission. In half of the IC 348 objects, 
the Ha profiles are broad and asymmetric, indicative of on-going accretion. Of 
these, IC348-355 (M8) is the lowest mass object to date to show accretion- like 
Ha. Three of our ~ M6 IC 348 targets with broad Ha also harbor broad 01 
(8446A) and Call (8662A) emission, and one shows broad Hel (6678A) emission; 
these features are usually seen in strongly accreting classical T Tauri stars. We 
find that in very low mass accretors, the Ha profile may be somewhat narrower 
than that in higher mass stars. We propose that low accretion rates combined 
with small infall velocities at very low masses can conspire to produce this effect. 
In the non-accretors in our sample. Ha emission is commensurate with, or higher 
than, saturated levels in field M dwarfs of similar spectral type. Our results con- 
stitute the most compelling evidence to date that young brown dwarfs undergo 
a T Tauri-like accretion phase similar to that in stars. This is consistent with 
a common origin for most low-mass stars, brown dwarfs and isolated planetary 
mass objects. 



Subject headings: stars: low mass, brown dwarfs - stars: pre-main-sequence - 
circumstellar matter - planetary systems - stars: formation - techniques: spec- 
troscopic 
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1. Introduction 

Brown dwarfs, which straddle the mass range between stars and planets, appear to be 
common both in the field and in star-forming regions. Their ubiquity makes the question 
of their origin an important one, both for our understanding of brown dwarfs themselves as 
well as for theories on the formation of stars and planets. 

In the standard framework, a low-mass star forms out of a collapsing cloud fragment, 
and goes through a "T Tauri phase" , during which it accretes material from a surrounding 
disk, before arriving on the main sequence. There is ample observational evidence now to 
support many key aspects of this picture for young solar-mass stars. One of the defining 
characteristics of the T Tauri phase is a rich array of emission lines in the spectra, of which 
Ha is usually the most prominent (Joy 1945). The broad, asymmetric Ha line profiles, seen 
in "classical" T Tauri stars, are now thought to arise in the infall region, as are many of the 
metaUic emission fines (MuzeroUe, Hartmann & Calvet 1998 and references therein). 

Whether the same scenario holds for objects at and below the sub-stellar limit is an 
open question. It has been suggested, most recently by Padoan & Nordlund (2003), that 
brown dwarfs form in the same way as more massive stars, via 'turbulent fragmentation'. 
Reipurth & Clarke (2001) proposed an alternate scenario, which has been further explored 
through numerical simulations by Bate et al. (2002; 2003): in their model, brown dwarfs 
(and presumably isolated planetary mass objects or "plancmos") arc stellar embryos, ejected 
from newborn multiple systems before they can accrete sufficient mass to eventually fuse 
hydrogen. In this scenario, a stellar embryo competes with its sibfings in order to accrete 
infalling matter, and the one that grows slowest is most likely to be ejected through dynamical 
interactions. A key prediction of Reipurth & Clarke (2001) is that "sub-stellar equivalents to 
the classical T Tauri stars should be rather short-fived" because the ejection process fimits 
the amount of gas brought along in a disk. 

Studies of young sub-stellar objects could provide valuable clues to distinguish between 
these formation mechanisms. Therefore, we have undertaken a multi-faceted study of very 
low mass (VLM) objects in star-forming regions and their immediate circumstellar environ- 
ment. One of our key goals is to investigate whether some or all young brown dwarfs undergo 
a T Tauri-like phase, and if so how long that phase lasts. 

This paper reports on a search for accretion signatures in high-resolution spectra of 14 
objects near and below the sub-stellar boundary in IC 348 (~ 320 pc) and Taurus (~ 150 
pc). Our targets consist of M5-M8 sources in IC 348 from Luhman (1999) and the recently 
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identified M7-M8 brown dwarfs in Taurus from Martm et al. (2001)^. Combined with the 
spectra of 15 additional sources in the Upper Scorpius and p Ophiuchus star-forming regions 
(Jayawardhana, Mohanty & Basri 2002), we now have a substantial sample to characterize 
disk accretion in young brown dwarfs. 

In a subsequent paper, we will use these optical spectra to derive temperatures, gravities 
and masses for the current sample, as we did for Upper Sco and p Oph (Mohanty et al. 2003a, 
2003b). For now, we assume an age of ~ 1-2 Myr for IC 348 (Luhman et al. 2000) and ~ 
3 Myr for Taurus (White & Ghez 2001). If we further adopt the Luhman (1999) spectral 
type-temperature scale for low-mass pre-main sequence objects, then theoretical evolutionary 
models (Baraffe et al. 1998) imply masses of ~ 0.1 - 0.02 Mq for our M5-M8 objects. Thus, 
our sample spans a mass range from the stellar/sub-stellar boundary well into the brown 
dwarf domain (and possibly even into the planetary mass regime; Mohanty et al. 2003b). 

2. Observations and Analysis 

We obtained optical spectra of the target sample using the High Resolution Echelle 
Spectrometer (HIRES; Vogt et al. 1994) on the Keck I telescope on 2001 October 27 UT 
(CFHT-BD-Tau-2, 3 and 4) and on 2002 October 31 UT (the rest, plus a second spectrum 
of CFHT-3). With the 1.15" slit, the two-pixel-binned spectral resolution is R ~ 33,000. 
The instrument yielded 15 spectral orders in the 6390 - 8700 A wavelength range, with 
gaps between the orders, providing a variety of features related to youth and accretion 
activity. For comparison to our targets and to derive v sin i, we used M dwarf and M giant 
spectroscopic standards observed with the same HIRES set up. The data were reduced in a 
standard manner using IDL routines, as described in Basri et al. (2000). 

We derived rotational velocities {v sin i) of the targets by cross-correlating with a 'spun- 
up' template of a slowly rotating standard. Multiple spectral orders (~ 6), selected on the 
basis of an absence of strong telluric features, strong gravity-sensitive features, and stellar 
emission lines, were used in the cross-correlation analysis. Following White & Basri (2003), 
we used a combination of giant and dwarf spectra for the template (see detailed discussion 
in Mohanty & Basri 2003 and Mohanty et al. 2003a). 



^The spectral classification of the Taurus targets by Martin ct al. (2001) has been modified by Bricefio 
et al. (2002). We adopt the latter authors' types, which are consistent with Luhman's (1999) typing scheme 
for the IC 348 targets. 
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3. Results and Discussion 

Table 1 lists the derived v sin i, Ho; equivalent widths, and Ha full widths at 10% of 
the peak flux for our sample. Figure 1 shows Ha line profiles in our sample. Of the ten IC 
348 objects, five show broad, asymmetric Ha lines; three of those five also exhibit strong 
emission in 01 (8446A) and Call (8662 A) (Fig. 2). Of the four Taurus targets, CFHT- 
BD-Tau-3 and CFHT-BD-Tau-4 show very strong but relatively narrow Ha. We have clear 
detections of Li (6708 A) absorption in IC348-256, 286 and 353, and marginal detections 
in CFHT-BD-Tau-2 and 3. The other spectra are too noisy in this region for unambiguous 
detection (not surprising for these very faint, red objects). 

3.1. Disk Accretion 

The shape and width of the Ha emission profile is commonly used to discriminate 
between accretors and non-accretors among T Tauri stars (TTS). Stars exhibiting broad, 
asymmetric Ha lines with equivalent width larger than 10 A are generally categorized as 
classical TTS (CTTS), although this threshold value varies with spectral type (e.g., Martin 
1998). Recently White & Basri (2003; hereafter WB) have suggested that a full- width > 270 
kms~^, at 10% of the peak emission, is a better empirical indicator of accretion, independent 
of spectral type. 

By the latter criterion, two of our targets are accretors: IC348-165 (M5.25) and IC348- 
205 (M6), with 10% full- width > 300 kms^^ Moreover, both objects show broad (FWHM > 
100 kms~ ) emission in 01 (8446A) and Call (8662A); IC348-165 also exhibits broad emission 
in Hel (6678A). In CTTS, the presence of a broad component in these metallic emission lines 
is associated with moderate to high accretion rates (A'f> 10^^ M0/yr), such as in DG Tau 
and DR Tau (Hessman & Gucnther 1997). WB detected these lines in ClDA-1 (M5.5) and 
GM Tau (M6.5), two Taurus objects at the sub-stellar boundary that also exhibit broad 
Ha profiles and veihng, and inferred accretion rates (~ 10~^-^ M0/yr) similar to CTTS. In 
fact, MuzeroUe, Hartmann, & Calvet (1998; hereafter MHC) found a correlation between 
the accretion rates and the fine fiuxes in Call, 01 and Hel for their sample of CTTS, and 
argued that the broad components of these lines originate in magnetospheric infall regions. 
MHC also found that the 01 (8446A) line appears in emission only for those CTTS with 
accretion rates higher than a certain threshold value (~ 10^^ Mo/yr). Thus, our detection of 
broad 01 and Call emission in 1C348-165 and 1C348-205 (as well as broad Hel in the former) 
is significant, and implies that some VLM objects are accreting in a fashion analogous to 
CTTS. 
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We also find broad 01 and Call emission in a third object, IC348-415 (M6.5) (Fig. 2). 
The results of MHC quoted above then indicate that accretion is very hkely in this object 
as well. It is thus noteworthy that its Ha 10% full -width (~ 213 kms~^) is somewhat below 
the accretion cutoff (270 kms~^) defined by WB. We suggest that the latter cutoff may need 
revision for VLM objects, for the following reason. The broad component of the Ho; line is 
believed to arise in a nearly free-falling fiow between the inner edge of the disk and the stellar 
surface (e.g., Hartmann, Hewett, & Calvet 1994); the lower limit on the line-broadening is 
then set by the free-fall velocity (v//). For ages 1-3 Myr and substcllar masses (M < 0.08 
Mq), evolutionary models (Baraffe et al. 1998) yield Vff ~ 75-125 kms~^, implying Ha 
full- widths of order 150-250 kms~^, significantly less than the WB cutoff of ~ 270 kms~^. 
MuzeroUe et al. (2001) argue that Stark broadening can increase the Ha width beyond that 
expected from infall alone. However, this effect does not appear important at rates below 
~ 10~^ M0/yr; thus in mildly accreting brown dwarfs, the Ha line-width is likely to be 
set by Vff alone. For instance, in V410 Anon 13 (a Taurus M6 object with 5xl0~^^ 
M0/yr), Muzerolle et al. (2000) adequately model Ha with no recourse to Stark broadening. 
These authors also show that accretion rates < 10^^ Mq/jy are not expected to produce 
any optical veiling arising from the accretion shock. The WB cutoff, however, is defined in 
terms of optical veiling: only objects with veiling are adopted as accretors, and it is these 
that have Ha full- widths > 270 kms"-*^. Our foregoing discussion suggests that VLM objects 
with very low accretion rates can fail the WB accretion test: their Ha profiles are unlikely to 
be Stark-broadened beyond the Vff full-widths, which are less than 270 kms"^; furthermore, 
they will also not exhibit any veiling indicative of accretion. A case in point is V410 Anon 13 
itself: its Ha 10% full-width is ~ 250 kms~^ (slightly below the WB limit), and it shows no 
veiling (and no Call or 01 emission either), yet is an accretor. Similarly IC348-415, which 
does show Call and 01 but has an Ha full-width somewhat below the WB limit, is very 
likely to be accreting. 

Consequently, we adopt a 10% full-width ~ 200 kms~^ as our accretion cutoff. Given 
the rotational velocities in our sample (all v sin i < 50kms~^), such a width is hard to explain 
by rotational broadening and/or flaring, but is in agreement with the expected Vff. In that 
case, at least two other targets in our sample - IC348-382 (M5.5) and IC348-355 (M8) - are 
accretors as well, albeit with low accretion rates. Both show broad asymmetric Ha profiles 
(10% full- width > 200kms~^) associated with infall, even though the 10% width is somewhat 
below the WB limit and no Call or 01 emission is noticeable. If so, IC348-355 is the latest 
spectral type (lowest mass) accretor known to date. 

Putting together now all the young objects near or below the substellar boundary (~ M5 
and later) with published high-resolution optical spectra, we have 4 objects in p Ophiuchus, 
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10 in IC 348, 14 in Taurus and 11 in Upper Scorpius^. Of these, optical spectral signatures 
of accretion are found in 1 object in p Oph (GY 5), 5 in IC 348 (discussed above), 3 in 
Taurus (CIDA-1, GM Tau, V410 Anon 13) and 1 in Upper Sco (USco 75 (M6), adopting an 
accretion cutoff of ~ 200 kms~^ in 10% width). The vast majority of p Oph VLM objects 
were inaccessible to our optical spectroscopy because of significant extinction, presumably 
due to circumstellar as well as interstellar material. Thus, our (small) p Oph sample is 
heavily biased against possible accretors, and should not be used to estimate the accreting 
fraction. Considering the other three clusters, which are much less affected by this bias, we 
find that ~50% of the VLM objects show disk accretion at an age < 2 Myr (IC 348), 20% 
at age ~ 3 Myr (Taurus), and ^10% by ~ 5 Myr (Upper Sco). While there are uncertainties 
in the cluster ages, IC 348 is likely to be younger than Taurus and Upper Sco. Thus, we 
appear to be seeing a decrease in the fraction of accreting young sub-stellar objects with 
increasing age. 

A decrease in brown dwarf disk frequency with age is also seen in measurements of 
infrared excess. While a large fraction -~60-80%- of brown dwarfs in the < 3 Myr-old 
p Oph, IC 348, and Trapezium clusters show near-infrared excess indicative of disks, the 
fraction appears to be lower in the somewhat older a Orionis (~5-8 Myr) and TW Hydrae 
(~10 Myr) groups (Jayawardhana et al. 2003; Jayawardhana, Ardila & Stelzer 2002; Muench 
et al. 2001, Liu ct al. 2003). The disk fractions as measured by near-infrared excess are 
roughly comparable between TTS and sub-stellar objects in the same cluster. The evidence 
to date is consistent with similar timescales for inner disk dissipation in young brown dwarfs 
and low- mass stars (Jayawardhana et al. 1999; 2001; Haisch, Lada & Lada 2001). 



3.2. Rotation and Activity 

In field ~ M5-M8 dwarfs, a saturation-type rotation-activity connection is seen: stars 
rotating above a threshold velocity (~ 5 kms^^) exhibit a saturated level of chromospheric 
Ha flux, while slower rotators show a range of (smaller) fluxes (Mohanty & Basri 2003). We 
will explore this issue for our young VLM objects in a future paper, after deriving effective 
temperatures so that observed Ha equivalent widths can be converted to fluxes. However, 
we note here that the Ha widths in our sample of non-accretors are at least as large as in 
saturated field dwarfs of similar spectral type (~ 3-loA). Indeed, in three cases - IC348-256 
and CFHT-BD-Tau-3 and 4 - the equivalent widths are many times higher than usually seen 



^Most from this paper and Jayawardhana, Mohanty & Basri (2002); 9 in Taurus from WB and 1 from 
Muzerolle et al. (2000). 
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in quiescent M dwarfs, and are reminiscent of intense chromospheric flaring. Interestingly, 
CFHT-BD-3 shows very strong Ha emission in two observations separated by a year. All 
these facts indicate that chromospheric activity levels in young VLM objects are at least as 
high as, and possibly significantly higher than, in older field dwarfs of similar spectral type. 

4. Concluding Remcirks 

Our high-resolution optical spectra provide the most compelling evidence yet for a T 
Tauri-like accretion phase in many sub-stellar objects, qualitatively similar to that of their 
stellar counterparts. The lifetimes of inner disks around brown dwarfs do not appear to be 
vastly different from those of T Tauri disks. Thus, the observational evidence to date favors 
a similar mechanism for the formation of stellar and sub-stellar objects. 

While it is too early to rule out the ejected embryo hypothesis for the origin of at least 
some brown dwarfs, our results do not support it: the naive expectation of signficantly shorter 
lifetimes for sub-stellar disks is not borne out by the data. However, current numerical sim- 
ulations do not have sufficient resolution to determine whether a sub-stellar object with a 
small, truncated disk should be an accretor at a few million years (Bate et al. 2003). Future 
simulations should clarify this issue. On the observational front, longer-wavelength mea- 
surements with the Space InfraRed Telescope Facility and/or the Stratospheric Observatory 
For Infrared Astronomy could constrain the sizes of brown dwarf disks and thus determine 
if they are truncated as expected in the ejection scenario. Other important insights could 
come from investigating even younger objects, namely embedded (proto-) brown dwarfs, 
with high-resolution near-infrared spectroscopy. 
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Hawaiians, and express our gratitude for permission to observe from its summit. We also 
thank the Keck Observatory staff for their outstanding assistance over the past several years. 
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Fig. 1. — Ho; line profiles of the target sample. Spectra shown have been smoothed by a 
3-pixel boxcar; continuum and full width at 10% of the peak levels are marked by dotted 
hues. Thick black hues indicate accretors with broad Ha as well as Call and 01 emission; 
grey indicates probable accretors, based on the Ha profile-shape and 10% full-width. Insets 
zoom in on objects with low peak-fiux and noisy continuua, to clearly show the Ha detection. 
For CFHT-3, two spectra are shown, separated by a year; note the similarly strong emission 
both times. For CFHT-4, note the change in Y-axis scale; the peak flux in this object is 
much higher than in any other target in our sample. 
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Fig. 2.— Line profiles of 01 (8446A) and Call (8662A) in three IC 348 objects. Arrows 
indicate line-center. These spectra have been smoothed by a 4-pixel boxcar. 
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^Spectral types for IC 348 sources are from Luhman et al. (2003; for 367, 382 and 415) and 
Luhman (1999; for others). Spectral type for CFHT-BD-Tau-1 is from Martm et al. (2001). 
Spectral types for remaining CFHT Taurus sources are from Briceno et al. (2002), who use 
the same classification scheme as Luhman (1999) and somewhat different from Martin et 
al. (2001). For CFHT-2 and 3, their spectral types are respectively 0.5 and 1.25 subclasses 
earlier than those of Martin et al. (2001). 

''PseucJo-equivalent width; the plethora of molecular lines make it impossible to determine 
the true continuum in these stars. 

■^For 01 and Call, large absorption troughs near the emission lines (see Fig. 2) make 
line- width determination difficult; our errors are ~ ±15%. 



